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This  is  the  2nd  quarterly  Technical  Progress  Report  summarizing  progress  on  the  Sea-optimized 
Quantum  Key  Exchange  (SeaQuaKE)  project,  which  is  led  by  Applied  Communications 
Sciences  under  the  ONR  Free  Space  Optical  Quantum  Key  Distribution  Special  Notice  (13-SN- 
0004  under  ONRBAA 13-001). 

In  this  technical  report,  we  describe  details  of  the  hyperentanglement  source  architecture  we 
propose  for  use  in  a  maritime  quantum  communications  system,  including  some  of  the 
wavelength-dependent  considerations  for  technologies  needed  to  construct  such  sources.  In 
addition,  we  discuss  initial  simulation  results  using  MODTRAN  to  model  atmospheric 
transmission  loss  as  a  function  of  wavelength  under  a  range  of  visibility  conditions  and  aerosol 
model  scenarios. 
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The  objective  of  the  ONR  SeaQuaKE  project  is  to  optimize  the  performance  of  free-space  optical 
(FSO)  quantum  key  distribution  (QKD)  operating  under  challenging  maritime  atmospheric 
conditions.  In  particular,  a  modeling  framework  will  be  developed  to  guide  optimization  of  the 
system  operating  wavelength  in  order  to  maximize  throughput  and/or  transmission  distance  over 
a  wide  range  of  atmospheric  conditions.  The  framework  will  consider  the  major  components  of 
the  quantum  communication  system  including  the  transmitter,  quantum  channel,  and  receiver 
elements.  Applied  Communication  Sciences  (ACS)  will  focus  its  efforts  on  the  transmitter  and 
receiver  elements,  while  Stevens  Institute  of  Technology  (SIT)  will  focus  their  effort  on  the  free- 
space  channel. 

The  focus  areas  over  the  last  quarter  include  (i)  further  development  of  the  hyperentanglement 
transmitter  architecture  and  its  wavelength  dependencies  and  (ii)  modeling  transmission  loss  of 
the  atmospheric  free-space  channel.  These  two  technical  areas  are  being  led  by  ACS  and  SIT, 
respectively.  Progress  towards  developing  models  for  these  elements  in  the  quantum 
communications  system  is  described  below. 
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3.1  Source  Architecture  and  Modeling 

The  transmitter  in  our  system  will  be  hyperentanglement-based  in  order  to  enable  greater  system 
throughput.  We  will  eonsider  polarization  and  time-bin  for  the  entanglement  degrees  of  freedom 
due  to  superior  channel  characteristics  for  these  modes  as  well  as  good  coupling  into  and  out  of 
single  mode  fibers.  The  focus  of  our  source  effort  is  to  understand  the  critical  components  and 
key  performance  parameters  that  vary  as  a  function  of  system  operating  wavelength.  We  focus 
our  efforts  on  5  infrared  wavelength  bands  with  the  best  transmissivity  through  the  atmosphere: 
0.8,  1.3,  1.5,  3.5,  and  9  pm. 

The  hyperentanglement  source  can  be  broken  down  into  three  smaller  subsystems,  including  the 
pump  source,  time-bin  multiplexer,  polarization  entanglement  &  pair  generation  elements.  The 
approach  we  are  taking  is  to  survey  the  current  body  of  literature  on  these  sub-systems  to 
understand  state-of-the-art  performance  capabilities,  trends  towards  improved  performance,  as 
well  as  fundamentals  limits  of  specific  enabling  technologies.  It  will  be  critically  important  to 
understand  how  each  of  these  factors  varies  as  a  function  of  source  operating  wavelength. 
Ultimately,  these  factors  will  translate  to  a  range  of  higher-level  source  metrics  to  be  considered 
in  our  model  including  entanglement  quality,  pair  probability,  pulse  rate,  mode  quality,  and 
signal-to-noise  ratios. 


3.2  Channel  Model  and  Validation 

The  channel  consists  of  the  free- space  link  separating  the  transmitter  and  receiver  elements  of  the 
quantum  communications  systems.  We  are  assuming  the  quantum  signal  may  propagate  through 
a  challenging  maritime  environment  with  widely  variable  atmospheric  conditions  including  haze, 
fog,  clouds,  and  rain,  as  well  as  scattering  and  turbulence  effects.  We  will  also  be  considering 
additional  sources  of  noise  that  can  corrupt  the  quantum  signal  recovered  at  the  receiver, 
including  noise  introduced  by  sunlight  as  well  as  blackbody  radiation. 

The  first  channel  impairment  we  have  started  to  consider  through  modeling  is  atmospheric 
propagation  loss.  SIT  has  set  up  an  initial  MODTRAN  simulation  for  atmospheric  transmission 
modeling  based  on  version  5.2.1  of  the  software.  The  atmospheric  profile  selected  in 
MODTRAN  defines  the  temperature  and  ratio  of  major  gases  and  particles  assumed  in  the 
simulation,  which  include,  for  example,  H2O,  O2,  CO,  CO2,  and  N2O.  The  U.S.  standard 
atmospheric  profile  was  used  in  the  simulations,  where  the  mixing  ratio  of  major  gases  is 
illustrated  below: 
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Figure  1.  Mixing  ratio  of  major  gases  for  U.S.  Standard  atmospheric  profile  in  MODTRAN. 

Up  to  5  different  aerosol  models  were  considered  in  the  atmospheric  transmission  simulations, 
including  (i)  Navy  Maritime,  (ii)  Maritime,  (iii)  Rural,  (iv)  Urban,  and  (v)  Troposphere.  As  an 
example,  the  particles  distribution  for  the  Navy  Aerosol  Model  is  provided  below: 


Red:  overall  size  distribution 
Green:  dust-like  aerosol 
(0.03  |im  mean  radius) 

Blue:  Marine  aerosol 
(0.03,0.24, 2  |dm  mean  radius) 


dNydR 


Particle  size  distribution 


R  lum] 


Figure  2.  Particle  size  distribution  for  Navy  Aerosol  Model  in  MODTRAN. 

MODTRAN  simulations  were  performed  over  a  1-10  |dm  wavelength  range  (with  0.1,  1,  5,  15 
cm'^  resolution)  assuming  a  30  km  free-space  link  distance  at  an  altitude  of  10  m.  The  effective 
visibility  was  varied  from  5  km  to  50  km  by  scaling  the  particle  size  distribution  appropriately. 
Additional  MODTRAN  simulation  parameters  included  assuming  a  relative  humidity  of  50% 
and  a  wind  speed=5  m/s. 
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4  Results  and  Discussion 


4.1  Source  Architecture  and  Modeling 

We  have  begun  our  source  model  evaluations  by  considering  the  architecture  given  below  in 
Figure  3.  The  hyperentangled  source  is  broken  down  into  several  sub-systems  including  the 
pulsed  optical  pump,  the  time-bin  interferometer,  and  the  containing  fiber  loop  which 
generates  the  single  photon  pairs.  This  architecture  is  based  on  fiber-based  components  and 
connections.  Optical  fiber  interconnects  are  relatively  stable  and  maintain  the  photon  pairs  in  a 
single  spatial  mode,  which  is  critical  for  ideal  propagation  through  the  atmosphere,  x*^^  materials 
are  chosen  initially  due  to  their  established  ability  to  be  integrated  with  single  mode  fiber  and  the 
ease  with  which  they  can  be  utilized  for  spontaneous  four-wave  mixing.  We  anticipate, 
however,  that  our  future  efforts  will  also  include  other  types  of  nonlinear  materials. 


5^(3) 

media 


Pulsed  interferometer 
laser  source 


Polarization 

interferometer 


Nondegenerate 
entangled  pairs 


Figure  3  Source  model  architecture  for  generating  hyperentangled  photon  pairs.  Sub-systems  will  be  chosen  and 
evaluated  in  terms  of  meeting  the  source  metrics  stated  in  the  text. 


Our  efforts  this  quarter  have  focused  on  the  identification  of  ideal  optical  pump  sources,  which 
will  ultimately  drive  the  nonlinearity  which  produces  the  single  photon  pairs.  We  have  identified 
several  critical  parameters  which  must  be  satisfied  by  the  pump  so  that  the  desired  quantum 
communication  rates  and  distances  can  be  realized.  These  are  summarized  below  in  Table  1. 


Table  1.  Pump  Parameters 


Parameter 

Value 

Pulse  Rate 

1  GHz 

Time-bin  Separation 

<500  ps 

Pulse  FWHM 

10-100  ps 

Peak  Power 

~1  W 

A  pulse  rate  of  1  GHz  is  chosen  because  we  anticipate  it  can  be  reasonably  met  with  currently 
available  technology,  though  depending  on  channel  characteristics  and  detector  efficiencies, 
higher  pulse  rates  may  have  to  be  used  to  boost  key  rates.  This  pulse  rate  then  dictates  the 
requirements  for  the  time-bin  separation  (<  Vi  pulse  period)  and  the  pulse  temporal  full-width 
half-maximum  (FWHM)  so  that  there  is  no  pulse  field  overlap  at  the  output  of  the  source. 
Finally,  the  pulse  source  must  be  capable  of  outputting  sufficient  peak  powers  to  drive  the 
nonlinearity  hard  enough  to  maintain  sufficient  pair  production  probabilities,  |l.  High  p  will 
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ensure  good  key  rates,  though  a  limit  of  approximately  0.01  is  typically  assumed  in 
entanglement-based  QKD  systems  in  order  to  prevent  photon  splitting  attacks.  In  1550  nm 

systems  employing  dispersion-shifted  fiber  it  is  known  that  jUo^iyPLY  and  yL  ~  0.1  is  a 
reasonable  estimate.  This  implies  that  a  peak  power  of  ~  1  W  should  be  adequate  for  producing 
ideal  photon  pair  rates.  This  power  level  is  used  as  a  starting  point  for  evaluating  the  various 
technological  options  currently  available  for  optical  pumps.  Finally,  our  initial  technology 
survey  is  focusing  on  self-pulsing  sources,  such  as  mode-locked  lasers  (MLL),  due  to  low 
complexity  and  high  peak  powers.  MLL  powers  are  sufficiently  high  that  filtering  to  appropriate 
pulse  widths  and  or  passive  temporal  multiplexing  to  GHz  rates  is  possible.  The  first  step  in 
establishing  the  system  model  is  ensuring  that  technologies  are  available  in  the  wavelength  bands 
of  interest  which  are  capable  of  meeting  the  metrics  outlined  in  Table  1.  Our  current  results  are 
summarized  below  in  Table  2.  Many  of  these  options  are  commercially  available,  while  others 
are  likely  to  require  additional  technology  development  efforts. 


Table  2.  Pump  Technologies 


Wavelength  (p,m) 

MLL  Options 

0.8 

Er  +  SHG,Ti:Sa 

1.3 

Pr,  OPO 

1.5 

Er:glass 

3.5 

HoPr:ZBLAN 

9 

OCL,  Er+DFG 

Table  2  acronyms:  OPO  =  optical  parametric  oscillator,  Er  =  erbium,  SHG  =  second  harmonic  generation,  Ti:Sa  = 
titanium  sapphire,  Pr  =  praseodymium,  Ho  =  holmium,  ZBLAN  =  a  family  of  fluoride  glasses,  QCL=quantum 
cascade  laser,  DFG^difference  frequency  generation 

4.2  Channel  Model  and  Validation 

The  results  of  the  MODTRAN  simulation  of  transmissivity  for  a  30  km  free-space  link  assuming 
a  Navy  Maritime  aerosol  model,  U.S.  standard  atmosphere,  along  with  additional  simulation 
parameters  defined  early  in  Section  3.2  is  shown  in  Figure  4.  Under  high  visibility  conditions 
(e.g.  23,  50  km),  reasonably  good  transmittance  can  be  seen  in  1.5-2  pm,  2-2.5  pm,  3.5-4  pm, 
and  8-10  pm  wavelength  regions.  However,  as  visibility  is  reduced  (e.g.  5,  10  km),  these 
preliminary  simulations  results  suggest  significantly  lower  loss  at  longer  wavelengths,  with 
almost  5  orders  of  magnitude  increase  in  the  calculated  transmittance  at  10  pm  vs.  1.5  pm  when 
visibility  is  5  km. 

In  order  to  explore  the  sensitivity  of  atmospheric  transmission  to  different  aerosol  distributions, 
MODTRAN  simulations  were  performed  using  a  range  of  aerosol  models  including  (i)  Navy 
Maritime,  (ii)  Maritime,  (iii)  Rural,  (iv)  Urban,  and  (v)  Troposphere.  Simulations  of  both  low 
visibility  (e.g.  5  km)  and  higher  visibility  conditions  (e.g.  23  km)  were  performed,  and  the  results 
are  shown  in  Figure  1.  As  in  the  Navy  Maritime  environment,  there  is  a  general  trend  towards 
reduced  transmissivity,  particularly  at  shorter  wavelengths  (e.g.  1  pm),  when  visibility  is 
reduced. 
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Figure  4.  MODTRAN  atmospheric  simulation  of  30  km  free-space  link  with  different  visibility  levels. 

(Navy  Maritime  aerosol  model,  U.S.  standard  atmosphere,  additional  simulation  parameters  defined  in  Section  3.2) 


- RURAL  23km 

- RURAL  5km 

NAVY  MARITIME  23km 

- NAVY  MARITIM  5km 

- MARITIME  23km 

- URBAN  5km 

- TROPOSPHERE  50km 


Beam  Transmission  Test  using  different  Aerosols 


Wavolongth  [micron] 


Figure  5.  MODTRAN  atmospheric  simulation  of  30  km  free-space  link  with  different  visibility  levels  and  different 
aerosol  models.  (U.S.  standard  atmosphere,  additional  simulation  parameters  defined  in  Section  3.2) 
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4.3  Deliverables/Milestones 


Date 

Deliverable/Milestone 

Status 

June  2014 

Progress  Report  No.  1:  Year  1,  Quarter 

August  2014 

Progress  Report  No.  2:  Year  1,  2^^  Quarter 

November  2014 

Progress  Report  No.  3:  Year  1,  3^^  Quarter 

February  2015 

Progress  Report  No.  4:  Year  1,  4*  Quarter 

May  2015 

Progress  Report  No.  5:  Year  2,  Quarter 

August  2015 

Progress  Report  No.  6:  Year  2,  2^^  Quarter 

November  2015 

Progress  Report  No.  7:  Year  2,  3"^^  Quarter 

February  2016 

Progress  Report  No.  8:  Year  2,  4*  Quarter 

May  2016 

Progress  Report  No.  9:  Year  3,  P^  Quarter 

August  2016 

Progress  Report  No.  10:  Year  3,  2”‘*  Quarter 

November  2016 

Progress  Report  No.  11:  Year  3,  3^^  Quarter 

February  2017 

Final  Report 
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5  Conclusions 


In  the  last  quarter  of  the  SeaQuaKE  project,  ACS  has  gone  into  further  design  of  the 
hyperentanglement  source  architecture  and  the  requirements  on  the  associated  technologies  that 
would  allow  it  to  operate  at  different  free-space  transmission  bands  from  0.8  |J,m  to  10  |J,m.  In 
addition,  SIT  has  initiated  MODTRAN  simulations  that  point  to  substantial  benefits  of  operating 
in  MWIR  and  LWIR  bands  under  low  visibility  conditions.  Finally,  we  attended  an  ONR 
program  review,  which  was  held  at  Sandia  on  July  22-23,  2014.  No  problems  are  currently 
anticipated. 
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